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We study the decoupling transition of flux lattices in a layered superconductors at which the 
Josephson coupling J is renormalized to zero. We identify the order parameter and related corre- 
lations; the latter are shown to decay as a power law in the decoupled phase. Within 2nd order 
renormalization group we find that the transition is always continuous, in contrast with results of 
the self consistent harmonic approximation. The critical temperature for weak J is ~ l/B, where 
B is the magnetic field, while for strong J it is ~ l/y/B and is strongly enhanced. We show that 
renormaliztion group can be used to evaluate the Josephson plasma frequency and find that for 
^\ • weak J it is ~ 1/BT 2 in the decoupled phase. 
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I. INTRODUCTION 



In many aspects the behavior of the high-temperature superconductors in the mixed state differs from that of 
conventional superconductors. A combination of elevated critical temperatures, high anisotropy and short supercon- 
ducting coherence lengths considerably enhance the role of fluctuations on the vortex interaction, which result in a 
£^ ' noticeable change in the nature of the mixed state. 

The phase diagram of these layered superconductor in a magnetic field B perpendicular to the layers has been 
studied in considerable detail^. A first order transition in Y Ba-iCu-iOi (YBCO) and in BiiSr-iCaCu^O^ (BSCCO) 
has been interpreted as a flux lattice melting. The data suggests that the first order line terminates at a multicritical 
point, which for BSCCO is at B « 300 - 10 3 G and at temperatures T « 40 - hQKM, while for YBCQi it is at 
Bo ss 2 — 10T and To « 60 — 80K, depending on disorder and Oxygen concentration. At low temperatures as B is 
i ^ i ■ increased a " second peak" transition is manifested as a sharp increase in magnetization. The second peak transition is 
at B rs Bq starting from the critical point; it is weakly temperature dependent and shifts to lower fields with increasing 
CN| ■ disorder 2 . More recent data on BSCCO shows^, however, that the second peak line connects smoothly with the first 
order line; the presence of a multicritical point depends then on possibility that an additional line joins in the high 
field regime. Josephson plasma resonance studies^ have shown a significant reduction of the Josephson coupling 
at the second peak transition. The combination of reduced Josephson coupling and enhanced pinning indicates an 
unusual phase transition. 

The flux lattice can undergo a transition which is unique to a layered superconductor, i.e. a decoupling 
r^j- transition 8,9,10,11 . In this transition the Josephson coupling between layers vanishes while the lattice is maintained by 
the magnetic coupling. The theory of Daemen et al^S employed the method of self consistent harmonic approximation 
(SCHA) to find the decoupling temperature Td(B). The SCHA leads to a conceptual difficulty since it predicts that 
the average (cosc/>), where <j> is the Josephson phase, vanishes at T > Tj, Koshelev has showni^ that (cos</>) is finite 
at all temperatures and in fact accounts for the experimentally observed Josephson plasma resonance at high temper- 
atures. Thus the decoupling transition, as found by SCHA, needs to be reinterpreted. The correct order parameter 
is in fact a non-local one (Eq. [2H1 below) and corresponds to the helicity modulus^ or to the critical current. 

In the present work we expand our earlier presentation^ and study the decoupling transition by 2nd order renor- 
malization group (RG). Section II defines the model and its effective Hamiltonian. In section III we identify the order 
parameter and related correlations and show that the latter decay either exponentially in the coupled phase or as a 
power law in the decoupled phase. We also show that RG can be used to evaluate (cos0) and hence the Josephson 
■ plasma frequency and predict its behavior in the decoupled phase (section V). 

Our RG analysis (section IV) shows that the decoupling transition is continuous in the whole parameter range , in 



contrast with the SCHA result that it is 1st order above some critical value of the Josephson coupling^. We trace 
the latter result to a deficiency of the SCHA which does not allow for a scale dependent effective mass (section VI). 

Our analysis assumes that point defects like vacancies and interstitials (V-I) of the flux lattice are not present. 
In general, however, V-I defects are a relevant perturbation at decouplingi£*i&. The effect is rather weak for actual 
system parameters, as discussed in section VII. 

The role of point disorder on the decoupling has been studied separately^ii and in more detail in a companion 
article 18 . This study shows that disorder induced decoupling leads to an apparent discontinuity in the tilt modulus 
which then leads to an enhanced critical current and reduced domain size. These results are in accord with the 
Josephson plasma resonance&l and other data on the second peak transition. Of further interest is the effect of 
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columnar defects on the decoupling transitioniS, an effect which can further identify this transition. 



II. THE MODEL 



Consider a layered superconductor with a magnetic field perpendicular to the layers. The flux lattice can be 
considered as point vortices in each superconducting layer which are stacked one on top of the other. Each point 
vortex, or a pancake vortex, represents a singularity of the superconducting order parameter, i.e. the order parameter's 
phase in a given layer changes by 27r around the pancake vortex. These pancake vortices are coupled by their magnetic 
field as well as by the Joscphson coupling between nearest layers. 

The basic model for studying layered superconductors is the Lawrence-Doniach 20 Hamiltonian in terms of super- 
conducting phases (p n ( r ) on the n-th layer and the vector potential A(r, z), A z (r,z) (vectors such as A and r are 
2-dimensional parallel to the layers): 



H LD = ^[d 2 rdz [VxA(r,z)f + ^ ( ^ Vft (r)-A(r, Z )) 2 i( 

ab y, 



z — nd) 



WE/ 



d rcos 
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nd 



J(n-l)d 



A z (t, z)dz 



(1) 



where \ a \, and £ are the penetration length and coherence length parallel to the layers, respectively, d is the spacing 
between layers and $o = hc/2e is the flux quantum. The first term in Eq. Q is the magnetic energy, the second is 
the supercurrent energy while the last one is the Josephson coupling where J is the Josephson coupling energy per 
area £ 2 between neighboring layers. The model Eq. qualifies as the standard model for layered superconductors. 

In this section we outline the derivation of an effective Hamiltonian, which is the basis for RG analysis in the next 
section. The partition sum for Eq. (|TJ involves integrating over both (f n (r) and A(r, z), subject to a gauge condition. 
Since A(r, z) is a Gaussian field (choosing the axial gauge A z (r, z) — 0) we can shift A — * A + 5 A where A(r, z) now 
satisfies the saddle point equations for its x, y components 



V x V x A(r, z) 



\ 2 

A ab 



E 



[f^Vp B (r)-A(r, -nd) 



(2) 



and then fluctuations in SA decouple from those of </>„(r). 
integration only on the 4> n {r) variables, 



The partition sum at temperature T involves therefore 



X><Mr)exp[-W IT] 



(3) 



with A(r, z) in Eq. given by the solution of Eq. for each configuration of (j> n ( r ) ■ Note that since Eq. is 
gauge invariant under A — > A — |^Vx(r, nd) and <y5 n (r) — > </2 n (r) — x(r, nd) one can in fact choose any gauge. 
We now decompose <p n (r) to 



Vn(r) = ^(i-)+E s «( r ') a ( r " r ') 
9 n (r) = <p° n (r) - ^(r) 



(4) 
(5) 



where <p n (r) is the nonsingular part of <£>n(r), a(r) = arctan(r2/ri) is the angle at r = (ri,^), s„(r) = 1 at pancake 
vortex sites and s„(r) = otherwise. The sum in Eq. Q is then a sum on r' being the possible vortex positions on 
the n-th layer, e.g. a grid with spacing £. 

Solving Eq. J5J) for A in terms of 8 n (r) and s n , substituting in Eq. yields after a straightforward analysis the 
equivalent HamiltonianpSS 



H 



(i) - 



TL V + T~Lj + Hf , 



(6) 



where TC V is the vortex- vortex interaction via the 3D magnetic field, Tij is interlayer Josephson coupling and Hf is 
an energy due to fluctuations of the nonsingular phase: 



(7) 



r,n r' ,n / 
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U J = [ rfMcos[^(r)+^( S „(r')-s„-i(r'))a(r-r')]-l) (M 

n *' r-' 



«/ = \j ^G- f \ q M6^M 2 (9) 

Here (q, k) is a 3D wave vector which is the Fourier transform to r, z with \q\ < l/£ and |fc| < 7r/d, and 

$2d 2 1 1 

Gv(q,k) = -SL^- (10) 

d sinhqd 
4A^g smh J ^ + sur ^ 
- , 167T 3 d 2 / 4A 2 b . 2 fcd\ 

<*/(«,*) = ^(l + ^sm^j. (12) 

Consider now a flux lattice with equilibrium positions R; (e.g. a hexagonal lattice) where I labels the flux lines. 
The actual positions of the pancake vortices deviate from the perfect lattice positions by u" on the n-th layer. The 
function s n (r) is then 

(r)= ( 1 if r = R, + up 
71 ' \ otherwise. 

The Fourier transform 

u(q, k) = 2^ u ™ cxp(iqR; + iknd) 

n , I 

identifies longitudinal u l (q, k) = q- u(q, k)/q and transverse u tr (q, k) = [q x z] ■ u(q, A;)/<? components of u(q, k) 

u(q, fe) = u l (q, k)q/q + u tr (q, k)[q x z\jq . 
where z is a unit vector in the z direction. The inverse transform is 

uj l = a 2 d/ u(q,fc)e- 4C l R! - lfcnd (13) 

where J qk = J BZ d 2 qdk / (2tt) 3 , J BZ is the q integration over the Brillouin zone (of area 4ir 2 /a 2 ) while |fc| < 7r/d, and 

a 2 is the area of the flux lattice's unit cell. Note that 9(q, k) involves much higher q components, up to the cutoff l/£. 

The next step is to consider small displacements in F v , Eq. JJJ; this is equivalent to assuming that temperature is 
well below the melting temperature T m . This expansion identifies the magnetic contribution to the elastic constants, 
leading to a Hamiltonian of the form 

=\{da 2 ) 2 [ [ q 2 c° m + ky 4i (k)]\u tr (%k)\ 2 

1 f 1 ^ d 2 ndk 
+H el {u l (q,k)} + - J ^_Gj\q 7 k)\9(q,k)\ 2 

-(j/e) E / d2r cos f^w + I>( r R < o - a ( r R u r +1 )i) 



(14) 



where W e z{w'(q, &)} involves elastic constants of the longitudinal modes and k 2 = -4 sin 2 t=. The elastic constants 
for the transverse modes, excluding Josephson coupling terms, are given by^i 

c ° 4(fc) = sd^fcf 111 i + +e41° (15) 

(16) 



" 66 ~ 16da 2 
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where Qq = 47r_B/$o = 47r/a i.e. ttQq is the area of a Brillouin zone, and 



(IT. 



If the Josephson term is expanded it would add more terms to the elastic constants, however, near decoupling the 
nonlinearity of the cos term is essential. 

We can estimate the condition T <C T m by evaluating the melting temperature via the Lindeman criterion ((w") 2 )o ~ 
c 2 L a 2 where cl is a Lindeman number of order 0.15 and the average is with respect to the elastic terms. The average 
is dominated by the softer transverse modes yielding T m h\(2ir 2 \ 2 ab / a 2 ) w c 2 r, i.e r is the temperature scale for 
melting. Melting in the absence of Josephson coupling was in fact studied^!, showing that T m is between r/8 and 
the two-dimensional melting temperature ~ 0.004r, approaching the latter at high fields a <C X a b- The condition 
Td -C T m is in fact satisfied only at a < A a & [Eqs. I|50lt)5l66[) below], hence we limit our solutions to < O.OIt. 

The next simplification, in accord with the condition of small displacements \u™\ -C a is an expansion of the relative 
phase in the Josephson term. 

5>(r R ( <) - a(r -R t - * - u?) • Va(r - R") = 6 n (r) (18) 

i i 

Since b n (r) contains all q in its Fourier transform it is useful to separate from it the components within the 1st 
Brilluin zone, 



b n (r) = b n {v)+B n (r) 

6 n (r) = 2md [ i3#L e -«l*-*»*( e ** - l) ^' fc) 

Jbz ( 27r ) q 



B n (r) = 2md V / ^ e -^(q + Q)T-^ (e ^ _ 1} u(q, k) ■ [* X (q + Q)] (lfl) 

The last term displays the Fourier transform of Va(r), which for Q = projects the transverse displacement. 
The lowest order effect of thermal fluctuations is obtained as an average with respect to the elastic terms, i.e. 

(cos[0 n (r) + 6 n (r)]) = exp[-^(r)) - \{bl (r)) ] (20) 

The only singular term in this average is due to b n (r) from Eq. I|19l) which in fact drives the decoupling transition, i.e. 
even if the fluctuations in uf are small, their effect on the Josephson phase can be divergent. The average (£? 2 (r))o 
is readily shown to yield w T/t, up to logarithmic terms, and since T <C T below melting its effect in Ea. (|20|l is 
negligible. 

We wish to integrate out also the high q modes of 0(q, k), i.e. momenta in the range Qq < q < l/£. The effect 
from Eq. I|20(l is a factor 



D 



- -hSl',!d 2 q dkG i (q,k)/(2nf _ 



which can also be neglected. 

The effective Hamiltonian is now defined with a momentum cutoff q < Qq and the corresponding smallest length 
scale is a. The Josephson term involves J^2 r — > (J/a 2 ) J d 2 r so that J is now the Josephson coupling energy per 
area a 2 , i.e. J = J(a/£) 2 . Since b n (r) depends only on the transverse modes all longitudinal terms are decoupled and 
we can consider an effective Hamiltonian for the transverse modes 

ni3) = \^(Gj\q,k)\9(q,k)\* + Gi\qMK<lM a ) 
q : fc 

- (J/a 2 )^ f d 2 rcos(6 n (r) + b n (r)) (21) 

where 

^ 1(9 ' fc) ^ 2 (2^ [c° 44 + q 2 c° 6 /k 2 z ] . (22) 
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Finally we shift 

<£n(r) = b n (r) + e n (r) 

so that #(q, k) can be integrated out leading to a Gaussian term i J2qki G b(<l, k) + Gf(q, fc)]~ 1 |^>(q, fc)| 2 . It is readily 
shown that Gf(q,k)/Gb{q,k) <C 1 for all q, k, except when both k < 1/A a 6 and q 3> k\ a b/a which has negligible 
effect since large k dominates the following integrals. Therefore we neglect Gf(q,k) and our final Hamiltonian is 

n = \Y. G b 1 ^ k )\M,k)\ 2 -(J/a 2 )Y J f d 2 rcos0„(r) (23) 
q,fc n •* 

Our task is then to evaluate the partition sum 

Z = [ V<t>e- HWIT . (24) 



III. GENERAL PROPERTIES 



In this section we identify the order parameter of the decoupling transition and related correlation functions. Also 
a scheme for evaluating the frequency of the Josephson plasma resonance is given. 

In the RG procedure, as detailed in section IV, a phase transition is found at a temperature Td such that at T > Td 
J is irrelevant (scales to zero under RG) while at T < Td it is relevant (scales to strong coupling). To identify the order 
parameter of this transition we consider first the Hamiltonian Eq. 1(23() in presence of an external vector potential 
in the z direction A n (r) which is z independent between layers. The relative superconducting phase on neighboring 
layers (f) n (r) is then shifted by (2ird/4>o)A n (r) and the Hamiltonian becomes 

n A = \j2 G b 1 ^ k M^k)\ 2 - (J/a 2 )J2 I d 2 rcos[^(r)-^A„(r)]. (25) 
q,fe n J ^° 

The Josephson current is a derivative of the free energy T — —T\nZ 

, . c dT 2irc J , . , , , . 2ird , . ... 
U r, n) = - 1 ^ TT -^ = j sin n r - — A n r )]) A . 26 

The linear response to A n (r) has a ~ cos(f> n (r) term as well as a nonlocal term from the expansion of Tiji in 

Jz (r,n) = - {^£\ ^[(cos0 n (r))A„(r) - ^ W d 2 r'(sin 0„(r) sin (r'))A„, (r')] (27) 

where averages (...) are in the A = system. The superconducting response is identified by a uniform A n (r) = A 
so that the superconducting phase ip n ( r ) acquires a uniform twist y n ( r ) — (2' K /(j)o)An. Note that the partition sum 
excludes y n (r) with a global twist, hence the A term cannot be transformed away. The superconducting response Q 
is defined by j n ( r ) = QA, i.e. 

- ^ J J[(cos0) - / d 2 ^(sin0o(O)sin0„(r))] (28) 



a 



where (cos0„(r)) is n and r independent is written for short as (cos0). This order parameter was identified by Li 
and Teiteli^ as a helicity modulus and by Koshelev 12 as the superconducting response (i.e. the zero frequency limit 
of his Eq. 3). This order parameter signifies breaking of gauge symmetry: at T > Td the irrelevancy of J implies 
diverging fluctuations in <p n (r), hence A n (v) has no effect on Z and Q = 0. At T < Td J flows to strong coupling so 
that 0n (r) has finite fluctuations near the energy minimum where </>n(r) = 0; this implies that a gauge transformation 
of A n (r) must be combined with a change in <fi n (r). The manifestation of this broken gauge symmetry is Q ^ or a 
finite Josephson current. 

The Josephson plasma resonance frequency oj p i is a significant probe&X of a possible decoupling transition. As shown 
by KosheleAsi^ the superconducting response at oj p i is dominated by the 1st term of Eq. 1)28(1 . so that qj p i ~ (cos</>). 
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We reconsider this relation in section V while here we present a scheme for evaluating (cos <fi) . We note first that a 
derivative of the free energy yields 

a 2 d 

- r W) (29) 



NL 2 dJ" 



where L 2 is a layer area and N is the number of layers. While the details of the RG procedure are not needed in 
this section, some general properties can be derived from the asymptotic RG transformation. A general RG changes 
the original lattice unit a to a renormalized one a R while J (a) — > J [a ) = J R ; additional parameters in J 7 (J) (not 
displayed here) may also be renormalized while T itself changes by 

<LF = -NL 2 f(a)da (30) 

In 1st order RG f(a) = (section IVA) but becomes finite in 2nd order (section IVB) so that /(a) ~ J 2 (a). Hence 
integrating Eq. H30JI yields a J dependent term for either a relevant or irrelevant y, hence (cos (ft) is finite in either the 
coupled or the decoupled phases. We consider in particular the decoupled phase T > Td and a system size i »a. At 
the final stage of RG when a R — ► L the renormalized J R becomes extremely small and one can safely use 1st order 
RG which has the form (section IVA) 

dJ/J = 2(1 - Z)da/a (31) 

where Z — Z(J,T) contains renormalizations due to higher order RG terms at shorter scales; note that Z — 1 at 
T = T d since for Z > 1 Eq. shows that J R ~ L 2 ^^ -> at L -> oo. E.g. if 1st order RG is used all the way 
from the initial values then Z = T/T^ where is Eq. (|44|l below. Integrating Eq. (|30|l yields 



F{J R ) - F{J) = -NL 2 I f(a')da' . 

J a 



(32) 



At the scale a R = L the system has just one degree of freedom so that the term ~ <? 2 </>(q, k) in Eq. iffify is absent, 
hence 

F(J R ) = -T\n[j** d4e jR ™*/ T ] N = -l T N[(^) 2 + O(^) 4 ] . (33) 

Thus F(J R ) - L^ 1 -^ at T > Td so that at L -> oo 

(cos0) =a 2 ^y f{a')da' + 0{L 2 - iZ ). (34) 

RG provides then an efficient method for evaluating lo p i via (|34() . 

We proceed now to study correlations. We note first the correction in Eq. (|34|l (cos^)^ — (cos</>) ~ L 2 ~ AZ (the 
upper limit in J L f(a')da! gives the same exponent as can be seen from section IV). This correction determines the 
decay rate of the correlation function, by considering the 2nd derivative d 2 J-/dJ 2 ~ (d/dJ) J T>(f> cos (f>o(0)e~ H jZ. 
This involves a cos</> ra (r) correlation as well as (cos0) 2 from dZ/dJ, hence 

d(cos(f>) 1 



dJ 



= Ta~2zZj d2 ^ cos MO) cos M*))c (35) 



where (cos <po(0) cos <p n (r)) c = (cos 4>o(Q) cos <pn( r )) ~ ( cos 4>) 2 vanishes at large r. To reproduce the finite size correction 
L 2 ~~ AZ at T > Td the correlation must decay as (cos </>o(0) cos 4>o(r)) c ~ l/r AZ . 

In the coupled phase a approaches the correlation length ^ at which J/T 1 becomes a strong coupling. The 
Hamiltonian can then be expanded as H ~ ^2 q (q 2 + £d~ 2 )l^( c l' ^)l 2 ' keeping just the q — > form. Hence for r > £<j the 
system becomes Gaussian and the correlations decay exponentially. We conclude that 

(cos0 o (O)cos0 o (r)) c ~ l/r 4Z T > T d 

- e" 1 "/^ T <T d . (36) 



This distinctive behavior of the correlations serves as an additional identification of the phase transition. 
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We finally examine the validity of the high temperature expansion which provides an easy estimate of (cos <j>) . We 
define (...)o as an average with respect to the J = system, so that 

\[M0)-M T )] )o= T —T- / o^„o ,u\ / JZ^2n 1 „2/^ ( 37 ) 



where the cqq term is absorbed into q 2 = 4 ln(a / dy/n) / X 2 b (at the dominant k — ir/d where c® 4 is significantly softer). 
The k integration yields the 1st order T° (Eq. 1441 below) so that 

([0„(O) - 0„(r)] 2 ) o = it [ 2n d " q 2/ 2 M ~ e-^ r ) = 8t\n(rq u ) r > l/q u 
J Trq 2 (l + q z /q£) 

= ~2tr 2 q 2 u \n{rq u ) r<l/q u (38) 

where t = T/T®. Using the Gaussian average 

(cos0„(O)cos0„(r)) o =exp{-([0„(O)-0„(r)] 2 ) o /2} (39) 
and expansion of Z to 1st order in J/T we obtain 

J f IT J 

(cos./.) = — J dV(cos0 n (O)cos0„(r)) o « Ta 2 q 2( 2 t - 1) ' ^ 4 °^ 

the contribution of r < l/q u which becomes comparable at large t is omitted. To check the validity of (|40|l we note 
that it can also be obtained from the exact relation IMol) taken as a perturbation in J, where (cos <j>) ~ J 2 is formally of 
higher order. However, this term is actually ~ J 2 L 2 N, hence the perturbation expansion formally fails at L, N — > oo. 
We note that within this naive perturbation expansion the correlation decays (incorrectly) to zero, 

< cos0"(r)cos0"(O) >~ (rq u )- 4t + 0[(J/T) 2 (rq u ) 2 - 4t ] (41) 

with a finite contribution to (cos0) in Eq. 14U|) . We note also that the next order terms in J, while decaying more 
slowly, are still convergent in Eq. (|40() for t > 1. In fact the result l|40(l . quiet remarkably, reproduces the 2nd order 
RG result (up to a numerical prefactor) as found below. The reason is that the decay of Ij41|l is actually correct for 
(cos </>o(0) cos 4> n { v ))c a s in Eq. (|36|l (with Z = t in weak coupling), which when substituted in (|35|l reproduces the 
form (|¥0|) . 

We summarize the salient features of the decoupling transition: (i) Relevancy of J corresponds to a broken gauge 
symmetry in the coupled phase, (ii) an order parameter that corresponds to (i) is Q of Eq. I|28|l or the Josephson 
current, and (iii) decay of correlations, such as (cos<^o(O)cos0„(r)) c , is either exponential in the coupled phase or a 
power law in the decoupled phase. 



IV. RG SOLUTION 



The decoupling transition is driven by the singular response of the Josephson phase 6„(r) in Eq. I|19|) . The 
singularity is due to the long range effect of an individual shift it™ which decays as Va(r — R") ~ [r — R™] -1 (see 
Eq. I18[l . The contributions of many such small displacements at a given location r result in a divergent response 
at the decoupling transition. The RG method is designed to handle such divergences and avoid the pitfalls of naive 
perturbation expansions. 



A. 1st order RG 



The RG method proceeds by integrating out slices of high momentum shells A — dA < q < A, and the momentum 
cutoff A is successively reduced; initially A = Q - The field ^(q, k) is then decomposed as <£(q, k) = C(q, k) + %(q, k) 
where £(q, k) carries momenta in the range A — dA < q < A while %(q, k) has momenta q < A — dA. In this subsection 
we analyze the phase transition by using the 1st order equation in J. Expansion of Eq. I|24|l to 1st order in J and 
averaging on the high momentum shell yields 

<cos(C„(r) + X r l (r)) > c = cos Xn(r)(l - \ < (C«(r)) 2 >) . (42) 
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Defining 



we obtain < (Cn(r)) 2 >= d J f^g(A,k)dA/A. Rescaling a 



q 2 T 
8irdg(q,k) 

— > a + da (da/a = dA/A) leads to 1st order RG, i.e. it 



identifies the change in the coefficient y of the cos term, with the initial value yo — J/T, as 

*=[1- 
V 



^9{A(x)M- 

ZTt X 



(43) 



where x = Qg/A 2 is in the range 1 < x < oo. For a continuous transition the limiting g(A — > 0, A;) can be taken and 
the vanishing of the right hand side in Eq. 143(1 identifies the decoupling transition temperature (which is independent 
of y ) 



T o _ 4a 4 



dk 

44 (ft) 



This defines the units of our temperature variable t, 



(44) 



(45) 



and Eq. (|4*5|l has then the form (f5Tf> with Z = t. 

The asymptotic solution of Eq. 1(43(1 is y(a;) = yox 1 ' 1 . For t > 1 the asymptotic y(x) vanishes on long scales x — > oo, 
hence the meaning of decoupling is that the Josephson coupling vanishes on long scales. For t < 1 the coupling y(x) 
increases, RG stops then when x reaches Xd = {£,d/a) 2 where strong coupling y(x) « 1 is achieved; this identifies a 
correlation length 



^«a(y ) 1/[2(t - 1)] 



(46) 



An explicit form for can be derived by noting the significant softening of (ft) at k > 1/a which implies that the ft 
integration in g(A, ft) is dominated by k w ir/d. Hence we replace in Eq. ((1511 ln[(l+fc 2 /<3g) / (l+£ 2 fc 2 )] — > 2 \n(a/dy/Tr), 
resulting in 



Gb(q(x),k) 
g(q(x),k) 



= da 



2 Axt sin 2 (fcd/2) 
2tsin 2 (fcd/2) 



so that 



and the parameter— <? are 



d | g 3 (A(x),ft) = r 



l/(4gx) 



(47) 
(48) 

(49) 



4-A 2 fc 



In 



d0r 



; ln(a/ds/i?) 



(50) 



Note that l/4gx = q 2 /q 2 L where q u is defined below Eq. ((37(1 . A continuous transition is determined by the x — > 00 
behavior, hence Eq. ((49(1 — > < and the critical point is at t = 1. We will examine below the possibility of a 1st order 
transition, hence in general we keep ~ 1/x terms. The solution for Eq. ((43(1 is then 



y(x) = y x 



l + 4ff 
1 + Agx 



(51) 



and the phase transition is indeed at t = 1. 
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B. 2nd order RG 



We proceed to study 2nd order RG. Our main objective in this subsection is to see if RG can reproduce a 1st order 
transition as proposed within the SCHA methodic. 

The RG procedure for a Hamiltonian of the type (|23|) has been derived in appendix A of Rcf 2 - up to 2nd order in 
J. In this process new terms in the effective Hamiltonian are generated so that the partition sum has the form 



2?0exp 



d 2 qdk , 



(2tt) 2 



d 2 r 



<V(q>*OI0(q>*Or + / ^{ycos[0„(r)]+vcos[<£„(r)+^„+i(r)]} 



where we define 



1 



87rd \g(q, k) 



+ ho + hi cos kd 



(52) 



(53) 



The new variables generated by RG, v(x), ho(x), h\(x), are cutoff dependent with initial values of «o(l) = ^o(l) = 
/ii(l) = 0. Note that the renormalization of h° and h 1 is equivalent to renormalization of c!j 4 . The recursion relations 
to 2nd order in y are2£ 

dy = [y(l-X )+yv(X +X 1 )]dx/x 
dh = [y 2 X Q + Av 2 (X + X 1 )]dx/x 

1 



[v(l - 2X - 2Xx) - -^X^dxjx 



dv 

dhi = [4u 2 (A + X 1 )]dx/x 
where X n , n = 0, 1, are ho, h\ and x dependent, 

dk cos(nfcd) 



X n = d 



2n l/g(x, k) + ho + h\coskd 



(54) 



(55) 



We have absorbed factors 7, 7' of order 1 in the definitions of y, v which depend on the cutoff smoothing procedure^; 
the initial value of y is then y(l) = 7'y. The equations l|54|l are to be integrated from their initial values. 
To first order in y we rederive Eq. I|43() above, 



dlriy = [1 — Xo(x, ho = h\— 0)]dx/x 



(56) 



Before presenting numerical solutions for Eqs. (|54(l . it is instructive to consider a simple situation where v and hi 
are neglected. This is a reasonable approximation since we eventually find that v < y and hi < Hq. The RG equations 
arc then 



dy = y(l - X (x))dx/x 
dho = y 2 X {x)dx/x 

where [using the approximate form for the In, as above Eq. I47j 

X (x) = d 



2lT l+l/(4gx) , , s 
2t sin^(fcd/2) + a 0\^) 



(57) 



(58) 



We consider first the asymptotic solution at x — > 00 in the regime where y is relevant, i.e. the coupled phase. 
Integrating the 1st equation of (|57|l we obtain 



ln^ 

yo 



1 f X Y(^ dx ' 

lux - / X (x ) — 
Ji x 



We claim that the 2nd term converges, hence the parameter s 



°°dx^=d 



dk 
2^ 



dx 



x + l/4g 



'{kIj2) +Xh °^ 



(59) 



(60) 
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defines the asymptotic form of y, 

V = yoxe~ s ■ (61) 

The RG equations are valid only up to y ~ 1, i.e. for a small yo up to a large Xd = (£d/a) 2 , but not infinite. Eq. 
(|60|l therefore assumes that by formally extending Eq. I|57|) to x — > oo the integration range between ^ d and oo is 
negligible. 

We note that the h$(x) term represents an additional mass term in the propagator Gp>(q, k), Eq. ll5.'.iH if h${x) ~ x. 
The xho(x) term in Eq. (|60[) is then ~ x 2 and the integral is convergent. In the SCHA method such a mass term 
serves as a variational parameter and an analogous equation to <|60[) is derived in the next section. The essential point 
here is that xh (x) is ~ x 2 only asymptotically while its scale dependence at finite x can be different, resulting in a 
different critical behavior. The scale dependence of the "mass" ho(x)/x is a feature which is beyond either 1st order 
RG or SCHA. 

To complete the argument we need to show that ho(x) increases with x justifying the convergence of l|60|) . Assuming 
that ho(x) increases with x, we can use Xq ~ 1/ho to yield from Eqs. (|57I61|1 that indeed ho ~ y ~ x. [For x < Xd 
more terms in the the denominator of Eq. I|58|l need to be kept, modifying the way ho increases with x.] 

We proceed now to solve the RG equations l|57|) in general form. Integrating Eq. I|58|l we obtain an explicit form 

*o = : i ■f(ho,x,t) (62) 

1 + tho + 

where 



4gx 



^• af '*) = H 1 -VTrlj 

and b = tho/ (1 + l/(4gx) + hot). The function f(ho, x, t) varies slowly between 0.84 < f(ho, x, t) < 1. 
The RG Eqs. (|57|) for y and ho become 

dy = (1 _ f{ho,x,t) to 

(1 + ^0+4^) X 

= y 2 f(h ,x,t) dx 
l + th + ^ x 

If the l/(4gx) term in Eq. (|63[) is neglected and f(ho,x,t) = 1 is taken, Eq. (|(j3|> is equivalent to the Kosterlitz- 
Thouless equations^ which show that the critical temperature is enhanced by a factor 1 + yo and that y ~ x when y 
is relevant. The presence of the 1/Agx term leads, however, to a more significant enhancement which is measured by 
a parameter p, 



/ Tyo [tyo~ , 

p= V^ = v¥ (64) 

measuring the ratio of tw o small parameters yo and g. It is also useful to express p in terms of the Josephson length 
\j = $ a\/rf/(4A ah \/j^), i.e. p = {Sn^^a/igXj). 

In appendix A we derive a form for the critical temperature which yields a proper p dependence for both small and 
large p, 

T d = (l + 2 P 2 )g T p«l (65) 
3 

T d = -^grp p > 1 (66) 

We have solved Eq. (|63J) numerically and the results for Td are shown in Fig. 1. The p-dependence is indeed 
significant, becoming linear at high p; in the range 1 < p < 1000 we find td = 0.95p . Thus Td ~ 1/B (Eq. I50fl at 
p C 1 while Td ~ (a/Xj)r ~ 1/yB at p ^> 1. We note also the significant enhancement in the value of Td in the 
latter case, i.e. T d > T°. We note that the RG expansion is valid for yo <C 1, which limits Eq. (|66|l (by inserting it 
in Eq. (JjU) to 1 < p < 1/g. 

At T < Td we find that the variable y first increases, then decreases, corresponding to weakening of the l/4g:z; term, 
and finally increases up to y w 1 at the scale x = Xd = (Cd/a) 2 where we should stop the renormalization process. 
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FIG. 1: The decoupling temperature td(p) = Td(p)/T® derived from the numerical solution of Eq. 1)63^ with 4g — 0.0001. 
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FIG. 2: The temperature dependence of the correlation length £d in terms of C = [1 + 21n(£d/a)] 1 near the decoupling 
transition temperature ta = Td/T® = 752.839. The numerical solution of Eq. 16311 uses 4g — 0.0001 and initial conditions 
yo = 0.01, ho = 0, . 

The temperature dependence of Xd is shown on Fig 2. We see that the correlation length £d — ► oo as T — ► which 
shows that the phase transition is a continuous one. We note that the early estimate of Glazman and Koshelev 9 of the 
decoupling temperature gives a result similar to Eq. (|66|l . They derive a condition for large fluctuations ((f) 2 ) which 
by itself does not prove a phase transition; furthermore the estimated critical temperature vanishes with p, i.e. it is 
incorrect at p <C 1. The large fluctuation condition is close in spirit to the SCHA method and is further discussed in 
section V. 

We present now numerical solutions of the full RG equations l|54Jl with Xq and X\ given by Ea. I|B5|I in Figs. 3, 4, for 
the same initial conditions as in Fig. 2 [yo = 0.01, ho(x = 1) = 0] and = 0, hi(x = 1) = 0. We choose Ag = 0.0001 
and these initial conditions since in this case the SCHA method (see section IV) yields a 1st order transition. 
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FIG. 3: The numerical solution of the RG Eqs. 1541 with 4g = 0.0001 and initial conditions yo = 0.01 and ho = hi — 0, 
projected on the y — ho plane for different temperatures near the temperature of the decoupling transition t d — T d /T d — 825.7. 
A higher curve corresponds to a higher t, i.e. the two lower curves have t <t d while the others have t > t d . 
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FIG. 4: Same as Fig. 3 projected on the y and v plane. 



Figures 3 and 4 show the projected solution on the y — h Q and y — v planes, respectively. We find that — 825.7; 
for t > td both the y and v variables vanish asymptotically, while at t < td both y and v are relevant. The resulting 
correlation length diverges at td, similar to Fig. 2, i.e. the transition is continuous. We have examined the transition 
also by varying the initial yo, vq; the correlation length was always found to diverge at T — * Td denning a continuous 
type phase transition. 

Another scenario for a 1st order transition is via changing the initial value of vq. Assuming a flow into strong 
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FIG. 5: The correlation length in terms of C = [1 + 2 ln(£d/a)] _1 for t = 0.5 and p = 1 as function of p v = y/ t\vo\/2g with 
Wo < 0. The cutoff on the RG integration of Eq. J2U is V / V 2 + v2 = lo3 - 



renormalized y r , v r the free energy is dominated by 

T r = ^[-y r cos(</>„) - v r cos(0„ + <p n +i)] (67) 

The minimum at 6 = 6 n = 9 n+ i changes from 9 — to cosf? = —y r /Av r at large a negative v r , v r < —y r /A. This 
hints at a 1st order transition at some initial negative v; this is not the decoupling transition, but rather a transition 
within the coupled phase . 

To emphasize the asymptotic forms we have integrated the RG equations up to \/y 2 + v 2 = 10 3 and show in 
Fig. 5 the corresponding correlation length. For these parameters, the asymptotic v becomes negative above p v — 
v / t\v \/2g = 0.5385 (with vq < 0). The curve in Fig. 5 is continuous at this p v since Xd is dominated by y 3> v in 
the asymptotic regime. One needs to further increase \vq\ until the asymptotic v and y values become comparable. 
In fact at p v = 0.65 we observe a marked change in slope for Xdipv) m Fig. 5. The asymptotic form y ~ x a is found 
with a decreasing from 1 as p v increases, saturating at a = 1/2 when p v > 0.65. Therefore, a 1st order transition is 
possible within the coupled phase, associated with the relative strength of the renormalized y and v variables. The 
transition occurs when the initial is sufficiently negative. 

V. PLASMA RESONANCE 

We show here the relation us p i ~ (cos</>) and then apply Eq. I|34(l to evaluate (cos</>). In presence of a weak 
time dependent electric field E(t) in the direction perpendicular to the layers the Josephson relation imposes a time 
dependent addition 8<j) n (t) to the Josephson phase. The kinetic energy has the form 

f E 2 3 egti 2 x ^ f 2 dd(j> n {t) 2 

n 

where eo is a dielectric constant. Expanding the Josephson coupling — (J /a 2 ) J2 n J d 2 r cos[</>„(r)+50 n (t)] yields to 2nd 
order i(J/a 2 )(cos</>)[<5</> n (i)] 2 (the 1st order term has (sm<j)) = 0). This form neglects the possible time dependence of 
the Josephson term, e.g. dynamics of pancake vortices which are assumed to be slow on the scale of uj p i. In particular 
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the "phase slip" frequency was showni^ to be much smaller than u> p i. The plasma frequency is then 

2 167re 2 dJ 

We proceed to evaluate (cos 0) from Eq. 134|) . Within 2nd order RG the contribution to the free energy 2 ^ has the 
form of Eq. |JSD| i.e. dT = -NL 2 f(x)dx with 

f{x) = j'T{y 2 {x)X (x) + 2v 2 (x)[X (x)+X 1 (x)}}/x 2 . (70) 

We are mainly interested in T > Td where v <§; y (see Fig. 4) so that we can use the truncated Eqs. (|57|l . We use 
here the solution IjAlf) for cither p <C 1 or p ^> 1, written as 

where Z — tj and in general Z = Z(t, J). The asymptotic form is ~ x 1 " 2 so that Z{t = td) = 1 identifies the 
transition temperature, e.g. in 1st order RG Z — t and (|71|) reduces to Eq. 1)510. Using Xq(x) = 4gtx/(l + igtx/Z) 
and |J7TJ, Eq. ® finally yields 

/ ^ J 2 t Z 

where tg = T/t <C 1 is assumed, corresponding to our requirement that T is well below melting. In particular for the 
case p -C 1, which is relevant for BSSCO^, we have Z — t and 

<«■*> = V^a^i (73) 

which can also be derived directly with 1)510: at high temperatures 2t ^> 1 we obtain (cos0) = ^'Jt/KT 2 . note that 
Eq. (|73|l is reproduced by the high temperature expansion Eq. I|4U|) up to a prefactor 47'. While the latter expansion 
is in general deficient, for evaluating (cose/)) it is reasonable, as discussed below Eq. i|41|) . 

Our main result for the decoupled phase in p < 1 systems is Eq. (|73H . In comparison, the melted phase where 
individual pancakes arc uncorrclatcd has a correlation length of ~ a, so that (cos 0o(O) cos 4> n (v)) c « e~ r l a is a 
plausible guess, hence from Eq. (|35f) we have (cos <fi) « J/2T; this form with a prefactor of order 1 was confirmed by 
simulations 12 . Hence in the decoupled phase (cos</>) is larger by a factor ~ t/T; furthermore, since J ~ 1/B (J/a 2 
is B independent) in the melted phase (cos0) ~ 1/TB while in the decoupled phase it is (cos0) ~ l/[BT(2t — 1)] or 
~ 1/BT 2 not too close to decoupling. Thus the temperature dependence of uj v i can distinguish between decoupled 
and liquid phases. 

VI. THE SCHA METHOD 

In this section we derive the decoupling transition within the variational SCHA method, reproducing the results of 
ref^fi. In particular, this method results in a 1st order transition at large p. We compare the method to that of 2nd 
order RG and show where the deficiency of SCHA originates. 

The SCHA proceeds by searching for the optimal Gaussian Hamiltonian of the form 

Wo = \ I E sr'fa. fc )<Mq> W(q. k ) (74) 

so that G s (q, k) is determined by minimization of the variational free energy 

F var = F + < H — H > 

where the averaging as well as Fq correspond to Hq. F var is then 

Fvar = I I {— \nG s (q, k) + [G^ 1 (q, k) — G^ 1 (q, k)]G s (q, k)} i- e ^hS q ,k G =(?. fe ) (75) 

* Jq,k a 
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where the lnG s term corresponds to Fq and we have used < 0*(q, fc)0(q, k) >o= TG s (q, k). The last term in 1751) is 
the Josephson term with, 

< cos0(r) > = e ~^ T S ^ Gs{qM) = exp(-s v ) . (76) 

This defines a parameter s v ; the renormalized Josephson coupling is then Je _s ". 
Minimizing Eq. 1(75(1 yields 

G- 1 =G fe - 1 + ^exp(- SlJ ) (77) 

Using the form l|47|) and the variable a; with d 2 q/ (2tt) 2 — —dx/(ax) 2 , Eqs. 1(76177(1 reduce to a self consistent equation 
for s v 

/dk f°° 1 
L dx i "xv~ -- (78) 

1/1 2isin a (fc<Z/2) T 

This equation has exactly the same structure as that of Eq. (|60fl within 2nd order RG if the asymptotic form of 
RG variable ho(x) — > xyoe~ s is used. However the detailed ho(x) behavior is significant and can affect the critical 
properties; in fact, y and ho(x) are non-monotonic. 

We can perform the fc integration in Eq. 1(78(1 . neglecting the / type function as in Eq. 1(62(1 (i.e. replacing 
0.84 < / < 1 by / = 1) leading to 

Agt / % arctan V$&\L if D > 



1 + Agx + (4 9 y P Wes» \ -^ fZgg^g ^ if D < ^ 

where D = 16g 2 (Ap 2 exp(— s v ) — 1) and p is defined in Eq. (|64fl . 

If the transition is continuous then s v diverges near Td so that the effective Josephson coupling ~ e~ Sv — > 0; Eq. 
(|75|l then yields ^ = 1. The RG result shows instead a weak p dependence even at small p as in Eq. (|B3|l and Fig. 1. 

At a 1st order transition s v is finite; anticipating a large p, 4p 2 exp(— s v ) ^> 1 but D w (8gp) 2 e _s " <C 1, Ea. l(7^|l 
can be written as: 

Sve -^l 2 ^l. (80) 
2p 

The product s„e~ s "/ 2 is bounded by 2/e at s„ = 2, hence as temperature approaches Td from below s v increases up 
to s v = 2 but then jumps to s v = oo in the decoupled phase, i.e. a 1st order transition. The critical temperature is 
then 

t d = —p 1 « p « 1/g (81) 
7re 

This result is similar to that from RG, Eq. except that the slope is somewhat different. The significant difference 
is that RG yields a continuous transition even at large p. 
At even larger p, where p> 1/g, D ^> 1, Eq. l(T§|) yields 

s v e- s " = . (82) 

As above, s^e -5 " is bounded by 1/e at s„ = 1, hence a 1st order transition at 

t d = -gp 2 P » 1/5 (83) 
e 

The results td = 1 for weak p and Eq. 183|) for strong p are the results given in RefiiS. The intermediate range 
Eq. H81fl i s n °t mentioned there, though the plotted decoupling fields Bd(T) in their Fig. 1 are consistent with 
Bd(T) — 1/T 2 as from Eq. J51J. Furthermore, Eq. (g3J yields T d = (4gp) 2 r/4e > r which is incompatible with the 
requirement that T d is well below T m . 

It is interesting to note that the early estimate of Glazman and Koshele\^ of the decoupling temperature gives a 
result similar to 1)81(1 or (|66|l . Within this estimate the cos <j> term in Eq. 1(23(1 is expanded and the condition of large 
fluctuations ((f) 2 ) ps 1 with (<j) 2 ) = T f k [G^ 1 (q,k) + J/a 2 ]^ 1 yields Td- (This condition indicates decoupling, though 
by itself does not prove a phase transition.) The result is then the same as Eqs. ((76177(1 with s„ ~ 1, therefore it 
yields indeed a result close to that of Eq. (|STJl . 

Our main result in this section is to show the formal similarity between SCHA and 2nd order RG as well as an 
important difference, i.e. the mass term which is generated by RG is scale dependent. Both methods show significant 
enhancement of Td with increasing Josephson coupling, however the transition remains continuous in the RG solution. 
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VII. CONCLUSIONS 



In recent experiments on BSCCC 2 '^ the phase diagram has shown a number of low temperature phases. Most of 
these transitions are disorder driven by either bulk pinning or by surface barriers. In particular the possibility that 
the second peak transition is a disorder driven decoupling has been suggestedi^it The significant reduction of the 
Josephson plasma resonance at the second peak supports a decoupling scenario 6 '", however a conclusive signature for 
decoupling has not been shown so far. 

The signature of decoupling is that translational order is maintained, though with softer tilt modulous, while 
superconducting order is lost, i.e. Q — (Eq. I28|) and the critical current in the c direction vanish at T > Td- An 
additional signature is the power law decay of the Josephson correlation at T > Td- We find that the decoupling 
transition temperature is Td ~ 1/B at low fileds (T° of Eq. IKHjl while it changes to Td — 0.95Tjp ~ 1/yB at higher 
fields (a A j < A 2 ) with significantly enhanced temperatures Td ^> T$. 

We have shown that RG can be used to evaluate (cos </>) and hence the Josephson plasma frequency. In particular 
for weak coupling J, as in BSCCO, we find that uj p i ~ 1/[BT(2T '/Tj — 1)], in contrast to a ~ 1/BT behavior in the 
melted phase. This temperature dependence can serve to identify a decoupled phase. 

In the present work we assume that VT defects are not generated. Hence superconductivity is lost only in the c 
direction (i.e. I c = 0) while 2-dimcnsional superconductivity is maintained parallel to the layers. Our neglect of VT 
defects is in fact not justified, since they are generated at a lower temperature in the J = systemi£*i&, i.e. at T$/8. 
The true transition is a 3-dimensional one in which both decoupling and the defect transition coalesce, similar to the 
B = scenario 2 ^. The actual transition temperature T c is between Td and the defect transition. It can be estimated 
by the temperature at which the correlation lengths of the defects £de/ and £<j become comparable. Thus e.g., if 
£,def > £,d, the Josephson coupling is renormalized to strong coupling before the cos term feels the VT defects. Since 
£def ~ aexp(E c /T) (where E c as 0.2t is the pancake vortex core energy; if local lattice relaxation is included 2 ^ then 
E c » 0.04t) £def is exponentially large, and T c is close to Td unless J is extremely small. 

We expect for systems like BSCCO or YBCO that decoupling affects mostly superconductivity in the c direction. 
The current-voltage relation parallel to the layers is expected then to be nonlinear—, except at very low currents where 
the few VT defects would eventually lead to a linear Ohmic behavior. Similarly, the power law for the Josephson 
correlation would eventually, beyond the VT spacing ^e/ decay exponentially. 

In conclusion we have studied the meaning and critical properties of the decoupling transition. On the theory side, 
in our view this is one of the few transitions of vortex matter which is fully understood. It remains to be seen if 
experiment can also provide clear realizations for this type of transition. 
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APPENDIX A: EXPANSIONS FOR T d 



We present in this appendix an analytic expansion for the decoupling temperature Td within the reduced set of Eq. 
(|63|l with f(ho,x,t) = 1. The results show a significant enhancement when the parameter p in Eq. I|64|l is large. 

We consider T > Td where y flows to zero and ho reaches a finite asymptotic value , since the integration of dho 
in Eq. (|63|l is convergent. We assume that the integration of the y equation is dominated by /ij 3 , and will examine 
below the validity of this assumption. 

Integrating y(x) in Eci. (|63[) with ho — > /j.q° yields 

1 + 4.9(1 + ^8°) 



v = VoX \ t+i^+W) ) ■ (A1) 

We now substitute this y{x) solution into the ho equation l|63|) and solve for ho, 

p*Z 2 ( Z(2Z-l) 4Z(Z-1) 



t A {Z- 1)(2Z- 1) V (l + igtx/Z) 2 ^- 1 ) (l + Agtx/Z) 22 - 1 

(t/p-i)(2t/ P -i y 

(1 + Agtx) 2t /P 



(A2) 
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where is & l so included in Z = l/(l/t + /ig°) and the important parameter p is defined in Eq. (|64|l . Since y ~ 
at large x, Z = 1 defines the transition point. 

We can now examine the consistency of replacing ho(x) by in the y equation. The condition that ho(x) 
approaches /ig° relatively fast is that the x dependent terms in Eq. I|A2() are small at x > l/Agt, i.e. Z is small (but 
Z > 1), hence p and are large. It seems plausible then that our approximation is valid at p 3> 1. Alternatively, if 
ho(x) <C 1 then it has anyway a weak effect on the RG, i.e. the present derivation is valid at p -C 1. 

Eq. 1|A2|) shows that ho(x) converges to if Z > 1. We now substitute x — > oo in 1|A2J) and obtain a self consistent 
equation for h^, which for the variable Z becomes a cubic equation 

2 + 3/< Z 2 + ? It n Z ~ l i2 ,n = z3 + a ^ 2 + ai Z + a o = (A3) 



J (p/ty + 2/t {p/ty + 2/t (p/ty + 2/t 

This cubic equation has solutions Z > 1 only if the condition £)(t) < is satisfied, where 

D (t) = - ^a aia 2 + ^a a\ - j4^ a i a 2 + \ a l ■ ( A4 ) 
Therefore, D(t) = correspond to Z = 1 and defines the temperature of the phase transition Td as given in Eqs. 
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